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We study the problem of constructing mutually unbiased bases in dimension six. 
This approach is based on an efficient numerical method designed to find solutions 
£NJ to the quantum state reconstruction problem in finite dimensions. Our technique 

suggests the existence of previously unknown symmetries in Karlsson's non-affine 

^ (2) 

family Kq which we confirm analytically. Also, we obtain strong evidence that 
C/3 no more than three mutually unbiased bases can be constructed from pairs which 

^| contain members of some non-affine families of complex Hadamard matrices. 
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3 I. INTRODUCTION 

7— 1 In a Hilbert space of dimension d we can construct maximal sets of d + 1 mutually 

> 

unbiased (MU) bases when d is prime or prime power. Otherwise, analytical [IH1] and 

numerical [5H9] efforts to construct maximal sets fail and it is suspected that they do not 

q{ exist. Any pair of MU bases {B\, B2} is unitary equivalent to the pair {I, H}, where I is 
O 

{Sj the identity matrix and H is a complex Hadamard matrix. Therefore, the existence of 
mutually unbiased bases is closely related to the existence of complex Hadamard matrices. 
Given a pair of MU bases {I, H} the problem to find the complete set of vectors MU to 
both of them can, in principle, be solved by considering the Buchberger's algorithm [10J. 
However, when H belongs to a special family of complex Hadamard matrices, known 
as non-affine, even 16 GB were not enough memory for the algorithm to terminate and 
identify the solutions [IT] . 
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This paper presents a new method that numerically solves this problem with the same 
computational cost for both affine and non-affme families. Our method is not an algorithm 
because it does not stop with a definite answer but it converges very quickly. 

This work is organized as follows: In Section II, we briefly introduce complex Hadamard 
matrices and mutually unbiased bases. In Section III, we present the method to find MU 
vectors and we discuss its convergence. We successfully test our approach in Section IV by 
searching the known 48 vectors MU to the identity and the Fourier matrix in dimension 
six, and we obtain known triplets containing the identity and the Di^a matrix. Section 
V contains our main results: we could not find triplets of MU bases from considering 
complex Hadamard matrix belonging to the non-affme families K$ and Kg as well as 
other families contained in them. This study have allowed us to find new symmetries for 
the family Kq 2 \ 



II. COMPLEX HADAMARD MATRICES AND MUTUALLY UNBIASED 

BASES 

This section contains the minimal information about complex Hadamard matrices and 
mutually unbiased bases required to make the paper self-contained; more details can be 
found in [B] and for example. Two orthonormal bases {|y?fc)} and (|0z)} defined on 
a d-dimensional Hilbert space are mutually unbiased (MU) if they satisfy the property 

\M<t>l)\ 2 = ], (1) 

for every k, I = 0, . . . , d — 1. Maximal sets of d + 1 MU bases have been found in every 
prime [12] and prime power (13] dimension. A lower bound on their members can be 
established in the general case of d = p^ 1 . . -p 7 ^ 1 , where p^ 1 < ■ ■ ■ < p r ": it is known 
how to construct at least p^ 1 + 1 MU bases [H]. Here, d = p^ 1 . . .p r ™ is the prime power 
decomposition of the number d. In the particular case of d = 6 the lower bound is three, 
and this is the maximal number of MU bases attained so far. 

A square matrix is a complex Hadamard matrix if it has unimodular entries and or- 
thonormal columns. Such matrices exist in every dimension and the Fourier matrices 
represent the simplest proof of their existence. Complex Hadamard matrices have been 
extensively studied in recent years and they are very hard to find when d>5 0. We say 
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that two complex Hadamard matrices Hi and H 2 are equivalent (Hi ~ H 2 ) if there exist 
unitary diagonal operators Di, D 2 and permutation operators Pi,P 2 , such that 

Hi = PiDiH 2 D 2 P 2 . (2) 

A complex Hadamard matrix may belong to a continuous set of inequivalent complex 
Hadamard matrices called a family. A family of complex Hadamard matrices H(x) is 
affine if it can be cast in the form 

H(x) = H(0) o Exp(ii2(x)), (3) 

where R(x) is a real matrix for all value of the real parameter ) and Exp 

is the entry-wise exponential function given by 

Ex.p(iR(x)) lm = exp (iR(x)i m ), (4) 

and the symbol o denotes the Hadamard product (A o5) im = Ai m Bi m . The number s 
of independent parameters corresponds to the dimension of the family. If a continuous 
family of inequivalent complex Hadamard matrices is not affine we say it is non- affine. 

(2) 

For example, the families stemming from the Fourier matrix in dimension six (F 6 ) and 
the Dita family (Dg) are affine families, whereas the Karlsson families and 
are non-affine. The notation here considered is consistent with the catalog of complex 
Hadamard matrices presented by Bruzda-Tadej-Zyczkowski [15J. In this notation, the 
upper index denotes the dimension of the family and the lower index the dimension of 
the space where it is defined. If a complex Hadamard matrix H belongs neither to an 
affine nor to a non-affine family we call it isolated. In other words, it is impossible to 
obtain a complex Hadamard matrix inequivalent to H from infinitesimal perturbations to 
its entries. A set of MU bases is inextensible if no further orthonormal basis MU to every 
base of the set exists. A well-known fact is that any set of d + 1 MU bases is inextensible, 
and it is conjectured that every triplet of MU bases in dimension six is inextensible. We 
mention that the complete set of inextensible MU bases in dimensions d < 5 has been 
found using Buchberger's algorithm [16]. Also, a characterization of triplets of MU bases 
in d = 6 have been given if the second MU basis belongs to an affine family of Hadamard 
matrices [8]. In a recent paper, it has been analytically proven that given any triplet of 
MU product bases in dimension six it is not possible to find even a single vector MU to 
the triplet [TF] . 
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III. MU VECTORS AS FIXED POINTS 
A. The physical imposition operator 

In this section, we present a method that allows the numerical construction of all of 
vectors MU to a given pair of MU bases {I, if}. The iterative method used here allows 
us to efficiently generate highly accurate approximations to the solutions of the defining 
set of equations. The desired states are attractive fixed points of the physical imposition 
operator [18] which has been used previously to find those quantum states known as Pauli 
partners [19J. The problem has a unique answer if the given probability distributions are 
informationally complete; otherwise a finite or infinitely many number of solutions may 
exist. The physical imposition operator is useful for the problem studied here because 
the search for MU bases is a particular case of the quantum state reconstruction problem, 
namely to determine the quantum state of a physical system compatible with probability 
distributions obtained from actual measurements. 

To illustrate the concept of the imposition operator, let us assume that {|v?fc)} and 
are the eigenvectors bases of two observables acting in C d say, A and B, respectively. 
We suppose that two probability distributions {pj^ } and {p\ B ^} have been obtained by 
measuring of the observables A and B, respectively. For simplicity, the distributions are 
assumed to be given exactly, which is only possible when the ensemble of quantum states 
is infinite. To reconstruct a pure state |$) 6 C d compatible with the measurements, we 
need to find all solutions {|$)} of the following set of coupled non-linear equations: 

\(^m\ 2 = A A) , (5) 
m®)\ 2 = p\ b \ (6) 

where k, I — 0, . . . , d — 1. In order to find a solution we perform the following steps: 

1. Choose a quantum state |\l/o) & t random, which will be called the seed. 

2. Decompose the state |\l/o) in the basis {\^>k)}, 

d-1 

|*o) =J] c *i^>- ( 7 ) 

fc=0 
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3. Modify the amplitudes of the expansion coefficients Ck in order to impose the infor- 
mation about A, 

c k -+ \ZpFttt < 8 > 

I k I 

In the last step, we have replaced the amplitudes of the coefficients {ck} compatible with 
those of the observable A; note that we did not modify the phase factors Ck/\ck\ because 
we can not draw any conditions about them from the data {p^}- The physical imposition 
operator implements the transformations just described in one operation, 

when |^o) happens to be orthogonal to the state \(fk), for any k — 0, . . . , d — 1, we define 

(<p k \%)/\& k \V )\^l. (10) 

This operator is non-linear and its action on every quantum state is well-defined. The 
action of this operator on a randomly chosen state |\&o) can be interpreted as incorporating 
what we learn about the unknown state when the observable A is measured. In other 
words, the initial state |^o) has no information about the quantum system considered 
while the state T a ^a) contains all the information we have acquired by measuring A 
in the unknown state. Note that T a ^a) is idempotent because applying it once exhausts 
the information available about A. 

Next, we proceed in a similar way with the second observable B in order to try to 
find a solution of the set of Eqs. 
with the observable B, 

Unfortunately, the state 

|*i)=T B)p{s) (T Ap{ A)|*o)), (12) 

generally does not contain the complete information about both A and B: some of the 
information about A is destroyed when T B p (B) is imposed, which is a consequence of the 
commutation rule [A, B] ^ 0. If A and B commute, it is trivial to find a solution to Eqs. 



5p), defining the physical imposition operator associated 
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( 5p ), namely 

S = {|$)G^/|$)=r Ap( A,|$ ),V|$o}G^} (13) 

= eH/\*) = T B , p(B) \%),\/\V ) eH}, (14) 

In general, the state has the complete information about B and only partial infor- 
mation about A, so the composite operator T B)P {b)T a ^a) is not idempotent. Therefore, 
we can iterate the procedure just described and analyze the convergence of the sequence 

l*«> = (T Bip(S )T A>p( A))"^o), n e N. (15) 
It has been proven [19] that every solution of the system of coupled equations (5|6) is an 



attractive fixed point of T b ^{b)T a ^a). Moreover, this property also holds for a general set 
of observables A, B,C, . . . and probability distributions p( A \p( B \p( c \ . . . The iterations 
are robust under adding redundant information, and the sequences converge if and only if 
the probability distributions are compatible, in the sense that the Heisenberg uncertainty 
principle is not violated. 

The problem of constructing MU vectors is a particular case of the quantum state 
reconstruction problem just described. Let A and B be two observables with a pair of 
MU eigenbases B A and Bb- A vector is MU to the pair {B A , Bb} if it has equally weighted 
probability distributions with respect to both observables, that is, 

P k A) = \ and Pi B) = k, (16) 

for every k, I — 0, . . . , d—1. Interestingly, when the eigenvectors bases are MU, every basin 
of attraction is found to be of the same size, verified numerically in every prime dimension 
2 < d < 37 [19J, as well as in every simulation reported below for d = 6. This property 
indicates that the efficiency of the algorithm is maximal when the eigenvector bases of the 
observables are MU, because the number of randomly chosen seed states needed to find 
all solutions is minimized. This observation conforms with the idea that the redundancy 
of information is minimal when the observables have MU eigenvector bases. 

B. Convergence 



In order to analyze the convergence of the sequence |^ n ) defined in Eq.(15) we need 



to define a metric for quantum states. We want to determine when a solution given by 
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our method is a solution of the coupled system of equations given by Eqs. ( f5p| ) . Let A 
be an observable having the eigenvectors base {\(fk)}k=o,...,d-i and let \<p) and \ip) be two 
arbitrary quantum states. The distance between the probability distributions associated 
with the observable A in the states \<p) and can be defined by means of Hellinger's 
metric [20], 

d-i 

\2 



(17) 



fc=0 



This metric compares two probability distributions of the eigenvalues of a single observable 
and it is important to realize that this is a metric for probability distributions, not for 
states. In the present context, we need to consider more than one observable and the 
corresponding probability distributions. Therefore, we introduce the Hellinger metric for 
m observables, the so-called distributional metric 



^jr D %{\ct>)M)), 



(18) 



where D A j(\(j)} } is the Hellinger distance of the observable A J , defined in Eq.(17). 

In our study of MU bases, we will always start from a pair of bases Ba = {| k = 
0, . . . , d — 1} and Bb = {\4>i), I = 0, . . . , d — 1}. Now, assuming that |$) G C d is a vector 
MU to these bases, the expression 



T>A,B(\*n),\®)) 



lD\(\* n ),\$))+±D B (\*n),m) 



1 /d-1 d-1 

\ Vd\^ ^ 




(19) 



tells us how close the state |\? n ) is to being MU to Ba and Bb- We will say that a sequence 
has converged when 

ZU, B (|tfn>,|$))<0.01, (20) 

which means that the absolute error of the amplitudes is less than 8 x 10 -4 on average. 
Numerical simulations suggest that the absolute error of every amplitude of a solution 
is very close to the averaged error just mentioned. Given that the desired solutions are 
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(stable) attractive fixed points, our approximations must be close to the exact solutions of 
the problem. In the next section, we test our method by constructing known sets of states 
MU to a number of pairs consisting of the identity and a complex Hadamard matrix of 
order six. 

IV. TESTING THE METHOD: TAO, FOURIER AND DITA MATRICES 

In this section, we apply the approach described above to four cases which have been 
studied before, reproducing successfully known results. We will (i) search for states si- 
multaneously MU to the identity matrix I and Tao's matrix the only known isolated 
Hadamard matrix of order six; (ii) we will derive the complete set of vectors MU to 
the pair {I, F$} with the numerical results being, in fact, so accurate that we are able 
to deduce an interesting analytic result about this set; (iii) we are able to confirm that 
there are no quadruples containing members of the Fourier family F e (a, b), and (iv) we 
will search for states MU to the standard basis and members of the one-parameter Di^a 
family (c). The results presented in this section are summarized in the first three 
rows of Table fl] in Sec. IVTl 

(i) Tao's matrix S$: The pair {I,S'g ' ) } cannot be complemented by six orthogonal 
vectors to form a triplet of MU bases [21] . We are able to confirm this result by unsuc- 
cessfully searching for a third basis by means of the imposition operator. We found 90 
vectors MU to the pair {I, S^} but a third MU base cannot be constructed from them. 

(ii) Fourier matrix _F 6 : It is impossible to construct four MU bases which contain the 
pair {I, Fq} jU [H], |22l 123] ■ More specifically, it is known that 48 vectors exist which are 
MU to this pair of bases, giving rise to 16 different ways to construct a triplet of MU 
bases. 

We have been able to unambiguously identify 48 vectors MU to {I, F 6 }, and they agree 
with the known solutions [H El [2IH23]- A careful analysis of the numerical expressions 
revealed that the components of 12 of the vectors can be expressed solely in terms of sixth 
roots of unity while the remaining 36 vectors also depend on Bjorck's number, 




(21) 
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which is unimodular, and occurs as a, a*, a 2 , (a 2 )*, where * means complex conjugation. 
The analytic expression for the number a has been found from the numerical results by 
imposing the unbiasedness of a solution to the pair {I, F 6 }. 

It turns out that the 48 vectors can be grouped into three sets, each corresponding to 
one orbit under the Weyl-Heisenberg group. To see this, let us first define the displacement 
operators D p by 

D p = r PlP2 X Pl Z P2 , (22) 

with p = (pi,p 2 ) £ Z 2 ,, where X and Z are the shift and phase operators, respectively, 
defined by their actions on the states of the canonical basis, 

X\(p k ) = \(pk+i), and Z\ip k ) = u k \ip k ). (23) 

Three vectors generating the mentioned orbits under Weyl-Heisenberg group are given 

by 

vi = -^(l,i,w 4 ,i,l,iw 4 ), (24) 

v 2 = -y=(l,-i,uj 2 ,-i,l,-iu} 2 ), (25) 

t> 3 = —=(l,ia,a 2 ,—ia 2 ,—a,—i), (26) 
V6 

with u = e 2 ™/ 6 ; we noted that V\ and v 2 are eigenvectors of and D^^^ respectively, 

for every (/, = 0, . . . ,d— 1. Consequently, the two orbits generated by V\ and v 2 each have 
six elements only consisting of the so-called Gaussian states [6] which, in fact, can be 
written as product vectors if we swap components 2 and 5 [24] . We understand the origin 
of the Gaussian vectors very well: the eigenvectors of the operators {Z?(i )0 ), -D(o,i)> -D(i,i)} 
- that is, the bases consisting of the eigenvectors of X, Z and XZ - form a triplet of MU 
bases in any dimension |25j. Every eigenvector of XZ in dimension six is a member of the 
orbit generated by Vi, whereas the orbit generated by v 2 is its complex conjugated orbit, 
in agreement with the fact that v 2 — v*. Given a MU vector of the pair {I, F$} it is well 
known that its complex conjugate is also a MU vector in any finite dimension. Moreover, 
this is also valid in infinite dimension and it is related to Perelomov's conjecture about 
the existence of Pauli partners [26J. The origin of the vector t> 3 , however, is not clear to 
us; it is not an eigenvector of any displacement operator, and it gives rise to an orbit with 
36 different states. 
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Finally, in condensed form, the 48 MU vectors can be written as 

{£>(m,o)^iW 6 (27) 
{D M v 2 )}, &6 (28) 
{£W^sWz 6 - (29) 

Both the first and second set define a circulant matrix each while the last set give rise to 
six circulant matrices. These observations generalize to other dimensions d as follows. 

PROPOSITION IV. 1 Let \<f>) be a state MU to the pair {I, F d }, where F d is the Fourier 
matrix defined in dimension d. Then, the set {F> p \<ft)} p€Z 2 defines an orbit of MU vectors 
which has d elements if \<p) is an eigenvector of any operator D p and it has d 2 elements if 
\<f>) is not an eigenvector of any D p . 

The proof of the proposition is trivial because the pair of bases of C d defined by 
{I, Fd} is invariant under the action of the displacement operators. Also, the eigenvectors 
of D(p tftv -) are shifted cyclically under the action of Z^'^v) for every // = 0, . . . , d — 1 
and v' ^ v [27]. 

(iii) Fourier family _F 6 (a, 6): we have attempted to extend pairs of the form 
{I, Fft 2 \a, b)} for 10 5 randomly chosen values of a and 6, taken from the entire parame- 
ter range. In each case we found a triplet of MU bases and we could not find a single 
additional vector MU to it. 

(iv) Dita family D§(c): The Di^a family D^\c) [28] is an affine one-parameter family 
of complex Hadamard matrices in dimension six which is closely related to the Fourier fam- 
ily. Bengtsson et al. [6 J found two triplets of MU bases which extend the pair {I, D^\o)}. 
To do so, they used a modified "24th-roots program" which lists all orthonormal bases 
whose vectors have 24th roots of unity as well as the number b 2 = (— 1 + 2i)/\/5 (and its 
complex conjugate) as components. 

In the numerical simulations realized with the physical imposition operator we have 
considered 10 4 random choices of the parameter c. We found ten triplets of MU bases 
containing the pair {I, D(0)}, but only two of them are inequivalent, in agreement with 
[B]. It is not difficult to find the analytic form of the triples by an educated guess. The 
exact value of the number 6 2 was found again by imposing the unbiasedness of a solution 
to the pair {I,D(0)}. 
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Explicitly, the inequivalent triplets we find are {I, -D(O), H m }, m = 1, 2, where 
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(30) 



(31) 



and uj = e 2ni/24 . The MU vectors are given by the columns of the matrices Hi and H 2 . 
Interestingly, both of them are equivalent to a member of the Fourier family. We have 
verified that these analytical expressions are indeed solutions of the problem. 

The two inequivalent triplets were found among the first three triplets obtained numer- 
ically, and no other inequivalent triplet was found in the next 100 runs of our program. 
This represents strong numerical evidence that no more than two inequivalent triplets 
exist which contain the pair {I, D(0)}. Moreover, both triplets occur with nearly equal 
frequency: we found {I, D(Q), Hi} 48 times while {I, D(Q), H 2 } occurred 52 times, an 
observation which can be explained if one assumes that the basin of attraction of every 
MU vector has the same size. This apparent symmetry has been noticed so far in each 
imposition-operator search for MU bases, whatever the dimension d [T8| IT9]. 

Triplets of MU bases containing the pair {I, D^\c)} have also been found for many 
other values of the parameter c; none of the resulting triplets seems to allow for even a 
single further MU vector. 
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V. KARLSON'S NON-AFFINE FAMILIES 



The most interesting property of the method defined in Section III is that its compu- 
tational costs do not increase when we consider non-affine families. This advantage can 
be used to analyze the construction of triplets from a pair of the form {I, H} where H 
belongs to a non-affine family of complex Hadamard matrices. Little seems to be known 
about extending such triplets, so that the method presented here is the first efficient way 
to study them over the entire parameter range of the families. We have performed com- 
putations for pairs of the form {I, K$ } and {I, K§ }, and it seems that they cannot be 
extended to four MU bases (up to a possible null measure set of the parameters). More- 
over, many pairs can be extended to a triplet only for a non-trivial subset of parameters. 

(2) 

Let us start with the family Kq . 

A. Karlsson's biparametric family 

Karlsson has found a two-parameter non-affine family of complex Hadamard matrices 
in dimension six [21], which contains the families D§~\ and two subfamilies of the 
Fourier family. The Dit-a family D$ (t) is equivalent to the four corners Kq 2 \±tt/2, ±7r/2) 

(2) (2) (2) (2) 

whereas K\ (x, 0) ~ F 6 (x,x) and K^ '(0, x) ~ (Fq (x,x)) 1 . Also, Matolcsi family 
determines one of the diagonals, that is, Kq 2 \x,x) ~ Mq~\x). All these subfamilies 
are explicitly obtained from fQ in Karlsson's paper [29]. Note that a subset of the 
Fourier family and its transpose define the horizontal and vertical axes, respectively, of 

(2) (2) 

the parameter space of Kq . Also, the Fourier matrix F 6 (0, 0) is equivalent to the center 
AT< 2) (0,0). 

(2} 

The results of our attempts to extend pairs of the form {I, K 6 (xi, £2)} to triplets are 
presented in Figjl} In this figure, a black dot at the point (0:1,22) means that a triplet 
has been found. Previously known result about affine families indicates that a triplet 
of MU bases can be obtained in the full range of the family [HI El EDI EI], as far as we 
know. The evidence presented in Fig. [T] shows that triplets of MU bases only exist for a 
subset of parameters if the family is non-affine. In these simulations we have considered 



convergence of the sequences according to the upper bound 0.01 established in Eq.(20) 



and we have also considered the bounds 0.03, 0.05 and 0.08. In all the cases we found the 
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same results, which evidence the stability of the solutions. 

Furthermore, Fig. [T] clearly suggests the existence of new symmetries. The Fourier 
matrix F§ and the Dita family (t) seem to be privileged in the problem of constructing 

(2) 

triplets from a pair of the form {I, Kg \X\, 22)}: both of them are centers of symmetries 
in Fig. [l Let us prove all symmetries existing in the family defined as: 

/ill 1 1 1 \ 

1—1 Z\ —Z\ Z\ —Z\ 

1 z 2 —fi —z 2 f 2 — h* —z 2 f{' 

1 ~Z 2 ~Zif 2 * Z\Z 2 j\ —Zlfi Z\Z 2 fz 
1 z 1 —f'3* ~ z 2fC —fl -2:2/2 
\ 1 -Z 2 -Zxfi Z 1 Z 2 f 3 -Ztf 2 * Z\Z%f\ ) 

where Z\ = e lXl and z 2 = e tX2 , —n/2 < Xi,x 2 < n/2 and the four functions 



kP(x u x 2 ) 



(32) 



fi 
h 

h 
U 



f{+X 1 ,+X 2 ), 

f(+xx,-x 2 ), 

f(-x h -x 2 ), 
f(-x u +x 2 ), 



(33) 



are defined in terms of a single function, namely 



e i(i-K*)/a f cos / *i ~ 



f(xi,x 2 ) 



Karlsson has shown that 



— 1 sin 




1 



1 + sin(xi) sin(x 2 ) 4 



(34) 



f(xi + it,x 2 ) = z 2 f(xi,-x 2 ), 
f(x 1 ,x 2 + n) = z 1 f(-x 1 ,x 2 ), 



(35) 
(36) 



meaning that 



and 



Kf\ Xl +n,x 2 ) = K^(xx,x 2 )P 3A P 56 , 



(2), 



K$\x u x 2 + tt) = P 36 Pa 5 K^ , (x 1 ,x 2 ), 



(2), 



(37) 



(3? 
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(2) 

FIG. 1: Triplets of MU bases from Kg . A dot at the point \X\,x 2 ) indicates that the pair 

(2) 

{I, Kq '{x\,x 2 )} can be extended to a triplet of MU bases. 

where P34 and P$q are permutations matrices. Consequently, one may restrict both pa- 
rameters X\ and x 2 to the interval [— tt/2, 7r/2]. 

Inspired by the symmetries of the graph shown in Fig. [T] and taking into account 



Eqs.(33) and Eq.(34), we notice that 

Kf\ Xl ,-x 2 ) = P 36 P 45 K^{x 1 ,x 2 ), (39) 
Kf\-x 1: x 2 ) = P m P^Kf\x 1: x 2 )- (40) 

using the symmetry f(xi,x 2 ) = f(x 2 ,Xi), we also obtain 

Ki 2 \x l ,x 2 ) = Kj i 2 \x 2 ,x 1 ). (41) 



Eqs.(39) to (41) reveal the symmetry apparent in Fig. [Tj and we consider it unlikely 

(2) 

that any further symmetries exist. The family Kj- (xi,x 2 ) with —tt/2 < x\,x 2 < n/2 
is divided into eight triangles of the same area, each of them containing one copy of the 
complete family. Therefore, it is sufficient to consider values in the triangle X\ 6 [0, tt/2], 
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x 2 < Xi, i.e. the shaded area in Fig. [2j In this figure, we show that the Matolcsi's 
family is located on the both diagonals of the square. As we will show later, we 
can construct triplets of MU bases from which means that it should appear in Fig. 

However, the set {Mg x} is of measure zero within the family Kg(xx,X2)', since the 
parameters (aq,^) are chosen at random, the probability of observe it vanishes. 




(2) 

FIG. 2: The Karlsson's family Kg (x-\_,x<i) is equivalent to the gray triangle; note that — 5 < 



Xj < 2 ' j ' ~ 1> ^ - 



B. Karlsson's tri-parametric family 



A tri-parametric non-affme family of complex Hadamard matrices has been recently 
found by Karlsson [32J, reading explicitly: 



where 



i*2 Z\ 7i2 

^ Z^ ^Z^BZ\ -^Z/^AZi j 



(42) 



1 1 
1 -1 



A 




B = -Fo - A 



(43) 
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and 

Zi= ( 

\Zi -z, 

The entries of A are given by 




(44) 



A n = -l + i^(cos(0)+e-^sin(0)), 
A 12 = _i + i^(-cos(6/) + e^sin(0)), 



(45) 
(46) 



and the entries of Z,- are 



4 



(47) 
(48) 
(49) 
(50) 



M A \M B (zi)), 
z\ = M A (zf), 
z\ = Mb(zI). 

Here, M. denotes the Mobius transformation, defined by 

a a/ \ oiz- f3 

M ( z ) = 7^ *> 

p*z — a* 

with a a = Af 2 , (3a = Af±, and «b = Bf 2 , /3b = B n , and 9, (ft, ip e [0, 7r). 

This family contains the non-affine family Kg and it also contains the complete set 
of the so-called ^-reducible matrices. In dimension six, a complex Hadamard matrix is 
-^-reducible if it contains nine 2x2 submatrices that are Hadamard matrices. Let us 



(51) 



analyze an interesting particular case. It follows from Eq.(45) and (46) that the subfamily 
Kq (0, 0, ip) do not depend on the parameter <p. In this case, the Mobius transformations 



in Eqs.(48-50) turn into the identity irrespective of the value of z. Therefore, we obtain 



an affine one-parameter family 



i^ 3) (O,0,^) = (P i6 F 2 ®F 3 )oExp(iR(^)), 



(52) 



(2\ 

which is contained in the Fourier family F 6 . Here, the matrix Ryip) is defined by 
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e e 

(c) Triplets for (9, <p, tt/2) (d) Triplets for Kf ] (6, <j), 3tt/2) 

(3) 

FIG. 3: Triplets in Karlsson's family K$ (#, 0, ip) for some fixed values of 





• • 


• • 






• • 


• • 






• • 


• • 










•/ 



where • means a null entry. The permutation matrix interchanges rows 4 and 6. This 
subfamily allows a triplet for any ip G [0, ir), because it is contained in the Fourier family 
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? 6 (a, b), which admit a triplet for any value of a and b. 



3(V 



to Fig. 3(d) 



(3) 

Numerical simulations considering the family K & are shown in Fig. 
The plots suggest the existence of specific reflection symmetries but we have not been 
able to confirm them analytically even in the simplest case ■0 = 0. 



VI. SUMMARY AND CONCLUSIONS 



We have presented an efficient numerical method to construct sets of mutually unbiased 
bases in finite dimension. The main advantage of our method appears when non-affine 
families of complex Hadamard matrices are considered, where the standard method to 
solve coupled polynomial equations (Buchberger's algorithm) often stalls due to excessive 
memory requirements. Our method numerically solves the problem for non-affine families 
with the same computational cost as for affine families. 

To test our approach we first used it to determine the well known 48 MU vectors to 
the pair {I, F 6 }. We have been able to prove that they determine three orbits under the 
Weyl-Heisenberg group. Also, considering the Dit-a matrix we found two triplets that 
agree with a result obtained previously [B|. 

Table [I] summarizes the results obtained in this paper, indicating the maximal set of 
MU bases that can be constructed from pairs of MU bases associated with various families 
in dimension six. In the cases of affine families and the non-affine family Bg we found 
triplets in the entire range of the parameters, whereas the isolated matrix does not 
allow a triplet. This property of has been previously found by Brierley and Weigert 

In the simulation realized for the non-affine family B$ we have considered 10.000 
random choices of its parameter, and we obtained a triplet in every case. This is the only 
non-affine family where we found a triplet for every value. The non-affine family 
defined in the range t G (vr/2, n] U (37r/2, 27r] does not allow a triplet in its full range. We 
have realized three simulations considering 1,000; 10,000 and 100,000 random choices of 
the parameter t and we obtain the same results. That is, a part of the family does 
not allow us to construct a triplet of MU bases (see Table [T]). 

(2) (3) 

In the cases of Karlsson's families and we have considered 2 millon and 8 
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Pair 


Kind 


Maximal set of MU bases 




Isolated 


2 




Affine 


3, Vcg [-±,±] 




Affine 


3, Vo, 6 G [0, 2vr) 




Non- affine 


3, Vs G [ 7r,arccos(- 1 + v ^)] U [arccos( - 1 + v ^), vr] 


{I, Ml \t)} 


Non- affine 


< 


3 if t G [0.5309vr, 0.9157tt] 
3 if t G [1.53127T, 1.9163tt] 
2 otherwise 


{I,Kj } 2 \x 1 ,x 2 )} 


Non- affine 


{ 


3 in black regions of Fig. [T] 
2 in white regions of Fig. [l] 




Non- affine 


1 


3 in black regions from Fig. 3(a)|to|3(d)| 
2 in white regions from Fig. 3(a)|tc|3(d)| 



TABLE I: Maximal set of MU bases for some families. The last four cases are new results. 

millon random choices, respectively, sampling the entire parameter ranges of both families. 
We have shown that these two families extend to triplets at most. The property that a 
triplet of MU bases can be found only for a reduced set of parameters of a family is a new 
result presented here by the first time. In addition, we identified new symmetries that 

(2) 

reduce the range of the parameters of the family fQ . 

Finally, in our investigation of more than ten million complex Hadamard matrices be- 
longing to non-affine families we could not find a single vector being MU to a triplet. This 
evidence supports the conjecture that no more than three MU bases can be constructed 
in dimension six. 
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